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Abstract 

Local search algorithms use the neighborhood relations among search states 
and often perform well for a variety of NP-hard combinatorial search problems. 
This paper shows how quantum computers can also use these neighborhood 
relations. An example of such a local quantum search is evaluated empirically 
for the satisfiability (SAT) problem and shown to be particularly effective for 
highly constrained instances. For problems with an intermediate number of 
constraints, it is somewhat less effective at exploiting problem structure than 
incremental quantum methods, in spite of the much smaller search space used 
by the local method. 



1 Introduction 

Combinatorial search problems are among the most difficult computational tasks: the 
time required to solve them often grows exponentially with the size of the problem [|l6| . 
Many such problems have a great deal of structure, allowing heuristic methods to 



greatly reduce the rate of exponential growth. Quantum computers || y, [TT|, [H, [13 
p8| offer a new possibility for utilizing this structure with quantum parallelism, i.e., 
the ability to operate simultaneously on many classical search states, and interference 
among different paths through the search space. 
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While several algorithms have been developed |38|, ||, [L8|, |21], |22|, the extent to 
which quantum searches can improve on heuristically guided classical methods for 
NP searches remains an open question. Even if quantum computers are not useful 
for all combinatorial search problems, they may still be useful for many instances 
encountered in practice. This is an important distinction since typical instances 
of search problems are often much easier to solve than is suggested by worst case 
analyses, though even typical costs often grow exponentially on classical machines. 
The study of the average or typical behavior of search heuristics relies primarily on 
empirical evaluation. This is because the complicated conditional dependencies in 
search choices made by the heuristic often preclude a simple theoretical analysis, 
although phenomenological theories can give an approximate description of some 



generic behaviors pO, |26|, |33 



The simplest algorithms are the unstructured methods, such as generate- and-test, 
which amount to a random search among the states or a systematic enumeration of 
them without any use of prior results to guide future choices. An analogous unstruc- 
tured quantum search can improve on this classical method ITS . 



One structured approach to search builds solutions incrementally from smaller 
parts. These methods exploit the fact that in many problems the small parts can be 
tested for consistency before they are expanded to larger ones. When a small part is 
found to be inconsistent, all possible extensions of it will also be inconsistent, allowing 
an early pruning of the search. In such cases, the search backtracks to a prior decision 
point to try a different incremental construction. Analogies with this approach form 
the basis of previously proposed structured quantum searches [^T], [2^]. 

A second approach takes advantage of the clustering of solutions found in many 
search problems. That is, instead of being randomly distributed throughout the search 
space, the states have a simple neighborhood relationship such that states with a few 
or no conflicts tend to be near other such states. This neighborhood relationship is 
used by repair or local searches. Starting from a random state, they repeatedly select 
from among the current state's neighbors one that reduces the number of conflicts 
with the constraints. Such searches can become stuck in local minima but are often 



very effective |32], p6| . The problem of local minima can be addressed by allowing 
occasional changes that increase the number of conflicts p5[ . By comparison with 
incremental methods, local searches operate in a smaller search space, i.e., they have 
no need to consider the many small parts that might be composed into a full solution. 
On the other hand, they also disregard any information from these small parts that 
can guide the search toward a solution. 

In this paper, we present a new quantum search algorithm based on the neighbor- 
hood structure of the satisfiability problem. Such structure forms the basis of many 
local classical methods but has not yet been applied to quantum computation. Specif- 
ically the following two sections describe the satisfiability problem used to illustrate 
our algorithm and the ingredients of quantum programs. We then describe a quan- 
tum algorithm analogous to classical local search methods and evaluate its behavior. 
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Finally some open issues are discussed, including a variety of ways this approach 
can be extended. The appendices derive the properties of the algorithm's quantum 
operation and discuss approaches to the classical simulation of the algorithm. 

2 The Satisfiability Problem 

NP search problems have exponentially many possible states and a procedure that 
quickly checks whether a given state is a solution Jjl|. Constraint satisfaction prob- 
lems (CSPs) [[30| are an important example. A CSP consists of n variables, V\, . . . , V n , 
and the requirement to assign a value to each variable to satisfy given constraints. 
Searches examine various assignments, which give values to some of the variables. 
Complete assignments have a value for every variable. 

One important CSP is the satisfiability problem (SAT), which consists of a prepo- 
sitional formula in n variables and the requirement to find a value (true or false) for 
each variable that makes the formula true. This problem has N = 2 n complete as- 
signments. For A;-SAT, the formula consists of a conjunction of clauses and each 
clause is a disjunction of k variables, any of which may be negated. For k > 3 these 
problems are NP-complete. An example of such a clause for k — 3, with the third 
variable negated, is V\ OR Vz OR V3, which is false for exactly one assignment for 
these variables: {V\ = false, V2 = false, V3 = true}. With the formula expressed as a 
conjunction of clauses, a solution must satisfy every clause. Thus we can view each 
clause as a constraint that adds one conflict to all complete assignments that include 
assignments to the variables in the clause that make it false. By selecting the clauses 
to be distinct, the number of clauses m corresponds to the number of constraints in 
the problem. 

The assignments for SAT can also be viewed as bit strings with the correspondence 
that the i th bit is or 1 according to whether V* is assigned the value false or true, 
respectively. In turn, these bit strings are the binary representation of integers, 
ranging from to 2 n — 1. For definiteness, we associate the value of variable Vi with 
the bit in the integer corresponding to 2* _1 , i.e., V\ corresponds to the least significant 
bit. Each clause in the a fc-SAT formula specifies a combination of k bits that cause 
a conflict in any assignment containing all of them. For example, the clause V\ OR 
V2 OR V3 introduces a conflict in all assignments for which bits 1 and 2 are and bit 

3 is 1, i.e., assignments whose integer representation ends with the bit string 100. 

For bit strings r and s, we define \s\ to be the number of 1-bits in s and r A s to be 
the bitwise AND operation on r and s. Thus |r A s\ counts the number of 1-bits both 
assignments have in common, and can also be viewed as a bitwise dot product of the 
two assignments considered as vectors of bits. We also use d(r, s) as the Hamming 
distance between r and s, i.e., the number of positions at which they have different 
values. These quantities are related by 

d(r, s) = \r\ + \s\ — 2|r A s\ (1) 



3 



An example with n = 2 and two constraints, each involving a single variable, 
is V\ 7^ 1 and V 2 ^ 1. This problem has a unique solution: Vj. = 0, V 2 = 0, an 
assignment with the bit representation 00. This problem is also an instance of 1-SAT 
with the propositional formula V\ AND Vi- The N = 2 n = 4 complete assignments 
have the bit representations 00, 01, 10, and 11. 

The search algorithm presented below uses the average number of conflicts in an 
assignment, c avg . Each of the m distinct clauses in a fc-SAT problem introduces a 
single unique conflict in all complete assignments whose values for the variables spec- 
ified in the clause match those of the clause. Since there are 2 n ~ k such assignments, 
the total number of conflicts in all assignments is m2 n ~ h . Thus the average number 
of conflicts in an assignment is 

c avg = m/2 k (2) 

In evaluating average behavior of search algorithms, a common class of problems 
used is the ensemble of random /c-SAT problems. This ensemble is specified by values 
for the number of variables n, the size of the clauses k and the number of distinct^ 
clauses m. 

In using random /c-SAT to test our algorithm, we restrict attention to those with 
at least one solution. Specifically, to generate "random soluble problems", a set of 
m distinct clauses is randomly selected from the \l)2 k possible clauses | 35| . The 



resulting SAT problem is then solved with a classical backtrack search to see if it has 
any solutions. Only those instances with solutions are retained. 

Soluble problems are very rare among randomly generated instances with many 
constraints. Thus to allow comparison over the full range of m/n we generate a slightly 
different ensemble, random problems with a prespecified solution. To generate these 
problems, a random assignment is first selected to be a solution. Then, m distinct 
clauses are selected from among those that do not conflict with that prespecified 
solution. For /c-SAT, this leaves 

- 1) (3) 

possible clauses to select among, because the prespecified solution precludes any clause 
that exactly matches the assignments in that solution. 

Compared to random soluble problems, using a prespecified solution increases the 
likelihood of problems with a larger number of solutions, resulting in slightly easier 
problems both classically and for the quantum search algorithm described in this 
paper. However, at the extremes of m = and m = m max , both methods give the 
same results. Furthermore, when m/n is fairly large, almost all soluble problems have 
only a single solution, resulting in very little difference between these two methods of 
generating soluble problems. 




lr This ensemble differs slightly from other studies where the clauses are not required to be distinct. 
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3 Quantum Computers 



In spite of difficult conceptual issues |15], ^TJ, the properties of quantum mechanics 
needed to describe the search algorithm presented in this paper are relatively simple. 
Specifically, these properties are the rich set of states available to a quantum computer 
and the operations that manipulate these states. 

The state of a classical computer can be described by a string of bits. For exam- 
ple, the complete assignment of an n-variable SAT problem considered at some point 
in a search is given by n bits. Quantum computers use physical devices whose full 
quantum state can be controlled. For example |D| , an atom in its ground state could 
represent a bit set to 0, and an excited state for 1. The atom can be switched be- 
tween these states, e.g., with lasers of appropriate frequencies. Significantly, the atom 
can also be placed in a uniquely quantum mechanical superposition of these values, 
which can be denoted as a vector (ipoyipi), with a component (called an amplitude) 
for each of the corresponding classical states for the system. Similarly, with n such 
devices, the state of the quantum machine is described by a vector with 2 n ampli- 
tudes (i/; , V'l; • • • > VV-i)- The amplitudes are complex numbers that have a physical 
interpretation: when the computer's state is measured, the superposition changes, 
or collapses, to become the classical state s with probability (V'sl 2 - Thus amplitudes 
satisfy the normalization condition J2 S \4>s\ 2 — 1- I n contrast to a classical machine 
which, at any given step of its program, has a definite value for each bit, the quantum 
machine exists in a general superposition of 2 n such states. 

The change in the superposition brought about by measurement also means the 
amplitude values are not explicitly available for use in a computation. For example, 
an algorithm that attempts to test the phase of various amplitudes to determine the 
next operation would destroy the superposition. This observation greatly limits the 
types of quantum algorithms. Thus using this rich set of states requires operations 
that rapidly manipulate amplitudes in a superposition without explicitly measuring 
them. Because quantum mechanics is linear and the normalization condition must 



always be satisfied, these operations are limited to unitary linear operators [21] . That 



is, a state vector ip = (ip , . . . , VV-i) can be changed to a new vector ip' related to the 
original one by a unitary transformation, i.e., ip' = Uip where U is a unitary matrixF] 
of dimension 2 n x 2 n . In particular, this requires that the operations be reversible: 
each output is the result of a single input. In spite of the exponential size of the 
matrix, in many cases the operation can be performed in a time that grows only as a 
polynomial in n by quantum computers [^, [21]. 

An important example of such operations are reversible classical programs. Specif- 
ically, let P be such a program. Then for each classical state s, i.e., a string of bit 
values, it produces an output s' = P[s], and each output is produced by only a single 

2 A complex matrix U is unitary when U^U = I, where U' is the transpose of U with all elements 
changed to their complex conjugates. Examples include permutations, rotations and multiplication 
by phases (complex numbers whose magnitude is one). 
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input. A simple example is a program operating with two bits that replaces the first 
value with the exclusive-or of both bits and leaves the second value unchanged, i.e., 
P[00] = 00, P[01] = 11, P[10] = 10 and P[ll] = 01. Denoting the two bits by 
the binary variables x and y, such a program performs the operation x <— x <S> y, 
where <S> is the exclusive-or operation. The values of the variables x and y after this 
operation uniquely determine the original values, hence the program is reversible. 
When used with a quantum superposition, such classical programs operate indepen- 
dently and simultaneously on each component, i.e., P[(/0o> "01) • • • j VV-i)] produces 
a new superposition (i/;' Q , . . . , V^-i) where ijJ 8 , = ip s with s' = P[s\. This quantum 
parallelism allows a machine with n bits to operate simultaneously with 2™ different 
classical states. Continuing with the 2-bit exclusive-or example, if we started with 
the superposition (l/y/2, 1/2, 1/2, 0) the program would produce (l/y/2, 0, 1/2, 1/2). 
Equivalently, this program corresponds to the permutation matrix 



U 



These reversible classical programs amount to a permutation of the states. How- 
ever, unitary operations can also mix the amplitudes in a state vector. An example 
for n = 1 is 

1 (\ -1 
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^TfU J (4) 

This converts (1, 0), which could correspond to an atom prepared in its ground state, 
to (1, l)/v2, i.e., an equal superposition of the two states. Because these matrices 
combine complex numbers, it is possible for the combination of amplitudes to cancel, 



leaving no final amplitude in some of the states. This capability for interference |p~5 | 
distinguishes quantum computers from probabilistic classical machines. 

Designing a quantum search algorithm requires combining these available quantum 
operations in a manner that concentrates amplitudes in solutions. Classical search 
techniques can be a useful guide for such designs. That is, a quantum search method 
could, at least approximately, map amplitude from one state to others that would 
be considered after it by the corresponding classical method. In effect, this allows 
examining, and using interference from, all possible choices the classical search could 
have made, rather than being restricted to a single series of choices. The details of the 
particular problem being solved could be introduced by adjustments to the phases of 
the amplitudes based on testing states for consistency, a rapid operation for NP prob- 
lems. This technique, used with previous unstructured [jn| and structured El], 
quantum searches, neatly separates the design of the method that mixes amplitudes 
from any consideration of the detailed nature of specific problems. In the rest of this 
paper, we show how analogies with local or repair style searches can be used as the 
basis for quantum search. 
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4 Neighborhood Relations and Quantum Search 



Classical repair style searches consist of randomly selecting a complete assignment 
as the initial state, then repeatedly moving to one of the current state's neighbors 
until a solution is found or the number of search steps reaches a prespecified bound. 
Typically, at each step the search selects a neighbor that reduces the number of 
conflicts, if any such neighbor exists. If the search terminates without finding a 
solution, it is repeated from a new random initial state. 

The two key ingredients of these local search algorithms are identifying the neigh- 
boring states (i.e., those with a single different assignment) and then selecting from 
among these neighbors the next one to use. The identification of neighbors doesn't 
depend on constraints of the particular problem being solved. On the other hand, 
the selection process uses information on the problem instance such as the number 
of conflicts in the current state compared to the conflicts of the neighbors. In some 
versions of repair searches, such as simulated annealing the selection process also 
depends on the number of steps that have been performed (e.g., through a changing 
temperature parameter). 

The quantum algorithm introduced here operates on superpositions of all complete 
assignments for the SAT problem. Furthermore, instead of selecting only a single new 
state, each step of the algorithm considers all possible new states. This parallelism 
offers the possibility of improved performance by arranging for interference among 
many different classical search paths. A successful algorithm will adjust amplitudes 
along these paths so that search paths leading to solutions tend to combine construc- 
tively while those leading to nonsolutions combine destructively. 

We construct quantum analogs of the two key ingredients of classical repair searches. 
First, the identification of neighbors corresponds to an operation on the superposi- 
tion that moves amplitudes preferentially from a state to its neighbors. Second, the 
selection among the neighbors is performed by adjustments to the amplitude phases, 
based on the number of conflicts in the states and their neighbors, i.e., the same 
information that is used in the classical methods. 

To describe the resulting algorithm more explicitly, let ip™) De the amplitude of 
the assignment s after completing step j of the algorithm. A single trial of the search 
algorithm consists of the following parts: 

1. initialize amplitude equally among all the assignments, i.e., = l/\/~N for 
all s. 

2. iterate: for step j from 1 to J, adjust phases and then multiply by the matrix 



U described in §fO], to give 

vv w = E u rsPs ^- l) (5) 



where p s is the phase adjustment for assignment s as described in §P~2 
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3. measure the final superposition 



The initialization can be performed rapidly by applying the matrix of Eq. (|j) 
separately to each of the n bits. 

After J steps, the measurement gives a single assignment, which is a solution with 
probability 

^scln = EW) (6) 

s 

I 2 

with the sum over solutions. Here p(s) = W * s the probability to obtain the 
assignment s with the measurement. On average, the algorithm will need to be 
repeated T = l/P so i n times to find a solution. As with classical repair searches, this 
algorithm is incomplete: it can find a solution if one exists but not prove no solutions 
exist. The search cost can be characterized by the number of steps required to find a 
solution on average, i.e., C = JT. 

Completing the description of the algorithm requires specifying the number of 
steps J. Because the measurement destroys the superposition, it is not possible to 
evaluate P S oin explicitly at each step and use this information to determine when 
to halt. Instead J must be specified a priori based on readily available problem 
characteristics, such as the number of variables n and clauses m. The value of J may 
be changed from one trial to the next, which allows some exploration for suitable 
values ||. 

In our case, because the algorithm emphasizes mapping amplitude among neigh- 
bors, one might expect that about n steps are required to ensure an opportunity to 
move significant amplitude to solutions (since a solution is at most n neighbors away 
from any given assignment). However, the experiments reported below show fairly 
large amplitudes with somewhat fewer steps. In addition, instead of continuing the 
algorithm to maximize the probability to have a solution, a lower average search cost 
is sometimes possible by stopping earlier |J, a simple strategy for improving prob- 
abilistic algorithms p9| . Determining the number of steps that minimizes average 



cost for fc-SAT problems with given n and m remains an open problem, but at worst 
one could try the algorithm for all values of J up to n, resulting in a linear increase 
in the overall search cost. More sophisticated methods for finding a suitable number 
of steps to take have been proposed for the unstructured search algorithm || and 
similar techniques may be useful for this structured search as well. 



4.1 Emphasizing Neighbors with Amplitude Mixing 

For SAT problems with n variables, each state has n neighbors. Thus the matrix 
M giving the most direct correspondence with the classical repair searches would 
have M rs = whenever the assignments r and s are not neighbors and M rs = 1/ \fri 
whenever r and s are neighbors, i.e., d(r, s) = 1. However, this matrix is not unitary 
and so does not correspond to a quantum operation. 
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Instead we consider a more general class of matrices whose elements depend only 
on the Hamming distance between assignments, i.e., U rs = Ud( r , s ) which can be made 
unitary by appropriate choices of the values Ud for d — 0, . . . , n. To give a close 
correspondence with the classical method, we use the unitary matrix of this form 
with the largest possible value of u\, which governs the mapping between states and 
their neighbors. 

As shown in Appendix |A|, such matrices can be written as U = WDW where, for 
assignments r and s, 

W rs = -L(-l)\ rAs \ (7) 
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and D is a diagonal matrix of phases (complex numbers with magnitude equal to 
1) depending only on the number of 1-bits in the assignments, i.e., D rr = tu. The 
matrix W is unitary and can be implemented rapidly on quantum computers using 
a recursive decomposition ^ [Tj|. The largest possible value of u% is obtained by 
selecting to be 1 for h < n/2 and —1 for h > n/2, as shown in Appendix [A]. These 
values define the mixing matrix U used in the search algorithm. 

For example, when n = 2 and the states are considered in binary order, i.e., 00, 
01, 10, and 11, 

(I 1 1 -1\ 
11-11 
1-111 
V-l 1 1 1 J 

For instance, the value mapping the state 00 to its neighbors, 01 and 10, is 1/2, 
corresponding to the second and third elements of the first column of this matrix. 

Fig. [I] shows relative values of the mixing matrix elements for larger cases with 
even n. This mixing matrix has a regular pattern of real-valued entries. Specifically, 
each column, specifying the mapping from a particular assignment s to all the others, 
has a positive value for assignments r at Hamming distance and 1 from s, negative 
values for d(r, s) = 2 and 3, and so on. More generally, Ud is negative when d mod 4 
is 2 or 3, and otherwise is positive. We make use of this pattern in designing the 
phase adjustments, described below. 

When n is odd, the pattern for the matrix values is different: Ud = for even 
values of d and the values for odd d alternate in sign. That is, Ud is positive when 
d = 1 (mod 4), negative when d = 3 (mod 4) and otherwise is 0. 

Unlike the local classical search, there is some mixing between assignments that 
are not neighbors. Generally, these contributions are relatively small, and become 
smaller as n increases, for most assignments. While these additional connections 
make it more difficult to understand the algorithm purely in analogy with classical 
repair style searches, the long range couplings provide a possible mechanism to move 
amplitude away from local minima, i.e., inconsistent states all of whose neighbors 
have more conflicts. 
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Figure 1: Values of Ud/u\ for n = 8 (solid) and n = 20 (dotted). The values of u\ for n 
equal to 8 and 20 are 0.27 and 0.18, respectively. 



4.2 Phase Choices 

The phase adjustment for states, i.e., the values of p s in Eq. @, is the portion of 
the algorithm using information about the particular SAT problem to be solved. The 
key to designing an effective algorithm is selecting phases that, at least on average, 
move significant amplitude toward solutions. As with classical heuristics, the phase 
selected for a given assignment must be a function only of readily available informa- 
tion about that assignment and not require specific knowledge of the problem as a 
whole, such as the location of the solutions, that can only be obtained by extensive 
search. Furthermore, as described in §|3|, the phase choices can not depend on explicit 
knowledge of the amplitude associated with the state. For example, the ability to set 
p s based on the sign of ip s could greatly simplify the design of the algorithm, but such 
an operation cannot be physically realized. 

Following the analogy with classical local searches, we make use of the same infor- 
mation such searches use, namely the number of conflicts in an assignment and how 
it compares to its neighbors. We found several methods that concentrate amplitude 
into solutions to some extent. In this section, we describe two such phase choices. 
The first makes few assumptions about the structure of the problem so applies fairly 
well to a wide range of problems. The second makes more specific assumptions about 
the problem structure and is thereby able to make more effective use of the pattern 
of values in U, but only for a limited set of highly constrained problems. 



4.2.1 A Simple Phase Choice 



the phase choice 
using p s = 1 when s is 



In the previous structured and unstructured algorithms [18|, £1 
consists of inverting the phase of all inconsistent states, i.e 
consistent, and otherwise —1. The more recent structured method |22j also inverted 
phases of some consistent states depending on the step j of the algorithm. In all 
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these cases, the mapping matrix tended to concentrate amplitude toward those states 
whose phases were unchanged. As described in Appendix |A|, the matrix U of § |4.1| 
is related to those used these previous algorithms. Thus we examine using phase 
inversion as a simple approach to moving amplitude toward desired states. 

In our case, the mapping emphasizes moving amplitude among neighbors. Since 
neighbors tend to have similar number of conflicts, we want to combine the ability of 
phase inversion to move amplitude toward those states whose phases are unchanged 
with an emphasis on neighbors to take advantage of the large value of Ux. A simple 
way to do this is to only invert the phases of states with at least a certain number 
of conflicts, and gradually reduce this threshold during the search. This allows for 
amplitude to move between neighbors while also biasing the system toward states 
with fewer and fewer conflicts. In effect this technique attempts to combine amplitude 
coherently for states with fewer and fewer conflicts while also giving relatively large 
transfers to these states through the large mapping between neighbor states produced 
by the matrix U. Although this argument ignores the, sometimes counterproductive, 
contribution from non-neighbor assignments, it does serve to motivate the phase 
selection method. 

Specifically we select an initial threshold c start for the number of conflicts and then 
at step j the phase for assignment s is 



The algorithm operates for at most c start + 1 steps (after this point, Eq. (||) gives the 
same phase to all states so there are no further changes to the relative amplitudes). 

The exact choice for c sta rt is not critical as long as there are relatively many states 
with about that number of conflicts. A simple approach is to start with the average 
number of conflicts in assignments, as given by Eq. (0). Thus for our algorithm we 
take c st art = m/2 k . Although not used in this paper, a possible alternative for a 
given problem, would be to count the number of conflicts in a set of randomly chosen 
assignments to estimate the distribution of conflicts among assignments. This distri- 
bution could in turn suggest a more suitable starting threshold, as well as different 
amounts to lower the threshold with each step of the algorithm. 

4.2.2 Phases Based on Neighborhood Structure 

Another approach to the phase choices is based on the pattern of the mixing matrix, 
illustrated in Fig. |l|. Specifically, if we could identify a series of disjoint sets of 
assignments in which successive sets have a large number of neighbors and the final 
set consists of solutions, then we could attempt to move amplitude from one set to 
the next based on the large mapping between neighbors. To use this technique we 
would also need to be able to initially concentrate amplitude in the first of these sets. 

Identifying such a sequence of sets requires fairly regular relationships among the 
number of conflicts of neighboring assignments. We restrict attention to the case 




1 if s has more than c sta rt — (j — 1) conflicts 
1 otherwise 



(9) 
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where n is even. A similar argument can give a phase selection strategy for odd n 



based on the different pattern of the U matrix described in §[41 . 

To illustrate this approach, we consider the simple case of soluble problems with 
the maximum possible number of constraints, i.e., m = m max from Eq. (|3|). In this 
case, each assignment with c conflicts has c neighbors with c — 1 conflicts and the 
remaining n — c neighbors have c + 1 conflicts. Let Abetter ( s ) be the number of 
neighbors of assignment s that have fewer conflicts than s does. Then for these 
problems, Abetter (s) = c where c is the number of conflicts in the assignment s. 

Let Sb be the set of assignments with b better neighbors, i.e., {s| Abetter (s) = b}. 
The set So contains just the single solution of the problem. The maximum constrained 
problems have (j^j assignments in Sb- Since this quantity is largest when b = n/2, 
the initial superposition with equal amplitude in all states will have the largest total 
contribution for the assignments in S n / 2 - Thus we will select phases in a way that 
attempts to move significant amplitude from the assignments in Sb to those in 5&_i, 
starting with b = n/2. If successful, this process will concentrate amplitude in the 
solution after n/2 + 1 steps. 

We can improve this process by first increasing the amplitude in the assignments 
in S n /2- This is readily done using the pattern of signs of the matrix U. Specifically, 
the initial superposition gives a positive amplitude to every assignment. We can 
arrange for the first step of the algorithm to give a sum of positive contributions to 
every state in S n / 2 by selecting the signs of the phases for the states to equal those of 
the corresponding matrix elements in U. That is, assignments at Hamming distance 
d from those in S n / 2 , namely those in S n / 2 ±d, have p s equal to the sign of u^. 

For subsequent steps j > 1, we attempt to move as much amplitude as possible 
from assignments in S n /2-Q-2) to those in S n /2-(j-i)- To do this we make use of 
the fact that the mapping between neighbors has the largest matrix elements in U 
and these elements, u\ are positive. Thus, just considering the contribution from 
neighbors, we can expect a large contribution from assignments in S n /2-(j-2) to those 
in S n /2-(j-i) by selecting p s for these assignments to be positive and setting the 
phase for all other assignments to —1. This argument ignores the contributions from 
assignments at larger Hamming distances. However, if a large portion of the amplitude 
at step j is concentrated in assignments in S , „y 2 -(j-2), the other contributions to 
assignments in S n /2-u-\) will be small both because the corresponding Ud values are 
small and the amplitudes in the other assignments are relatively small. 

This discussion motivates selecting the phase for an assignment s as follows. For 
the first step, 

_ f -1 if |A^ start - Af bette r(s)| mod 4 = 2 or 3 , . 

^ s 1 1 otherwise 

where A^ start = \n/2\. For subsequent steps, i.e., j > 1, 

_ f 1 if Af start - Af better (s) = j — 1 or j — 2 
^ s 1—1 otherwise 



12 



with the algorithm terminating after at most j = N sta , r t + 1 steps. 

As suggested by the motivating argument, this choice of phases requires a close 
relation between the number of conflicts in an assignment and its neighborhood rela- 
tions. While such a relation holds exactly for the maximum constrained problem, it 
is only an approximation for problems with fewer constraints. For example, problems 
can have local minima, i.e., assignments for which Abetter = but which nevertheless 
have some conflicts and thus are not solutions. More generally, the neighbors of a 
given assignment can have various numbers of conflicts making it more difficult to 
move amplitude one step at a time to sets of assignments with successively fewer con- 
flicts. These characteristics introduce variations in the amplitudes, in effect causing 
some randomization of the phases. We observed the validity of the approximation, 
and its usefulness as a search strategy, was lowest for problems around the transition 
from over- to underconstrained problems near m/n = 4.2. Nevertheless, for highly 
constrained problems the behavior is close enough to that of maximum constrained 
cases to give the step by step shift of a large portion of the amplitude toward assign- 
ments with fewer conflicts. 

5 Quantum Search Behavior 

The behavior of this search algorithm was examined through a classical simulation. 
While these results are limited to small problems, they nevertheless give an indication 
of how this algorithm can dramatically focus amplitude into solutions. As a check on 
the numerical errors, the norm of the state vector remained within 10~ 10 of 1. 

5.1 Small Example 

We first illustrate the steps in this algorithm, for the trivial problem with two variables 
described in §[2|. For this problem the characteristics of the assignments and the 
amplitudes using Eq. @ are: 



assignment s 


00 


01 


10 


11 


number of conflicts 





1 


1 


2 




0.5 


0.5 


0.5 


0.5 




0.5 


0.5 


0.5 


-0.5 


1>P 


1 












with C s t ar t = 2/2 1 = 1. 

Initially the states are prepared with equal amplitudes. At the first step Eq. (^) 
inverts the amplitude of all states with more than one conflict, i.e., the state 11 in this 
case. Multiplication by the mixing matrix of Eq. (|8|) then gives the final amplitudes. 
In this case, after one step, all the amplitude is in the single solution state so P so i n = 1 
and the search cost is C — 1 step. This compares with an expected cost of 4 steps 
from random selection. 
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The neighborhood based phases choices of Eq. fllO| ) are based on each assignment's 
number of neighbors with fewer conflicts, and iV start = 1. For this problem, |iV start — 
-Wbetter(s)| < 1 so Eq. (|10a| ) gives no phase changes for the first step, i.e., ijjf' = 
Thus, for this small problem the attempt to concentrate additional amplitude in 
assignments with Abetter (s) = A^art ^ as no effect- For the next step, j = 2, Eq. (|10b| ) 
inverts the phase for the assignment with Abetter = 2, giving 



assignment s 


00 


01 


10 


11 


Abetter (s) 





1 


1 


2 


A'start — A/better (s) 


1 








-1 




0.5 


0.5 


0.5 


0.5 




0.5 


0.5 


0.5 


-0.5 




1 












Thus this method also gives P so i n = 1 but only after 2 steps, so the search cost is 
C = 2. 



5.2 Maximum Constrained Problems 

A particularly simple case is given by problems with the maximum possible number 
of constraints that still allows for a solution, i.e., m = m max as given in Eq. (|3]). 

In this case, classical repair style search methods can find a solution rapidly start- 
ing from any initial state, typically by removing one conflict at a time. Similarly, 
incremental classical methods will encounter conflicts immediately upon choosing 
any variable value not in the solution, thus allowing the solution to be found in 0(n) 
steps. 

As shown in Appendix p.2| , the classical simulation of the quantum algorithm for 



these problems runs rapidly. Thus we are able to explore the behavior of the algorithm 
with considerably more variables than is possible when there are fewer constraints. 

The scaling of the search cost for this extreme problem is shown in Fig. [2] on a 
log- log plot where a straight line corresponds to a polynomial growth in the cost. 
The expected search cost grows quite slowly and is approximately proportional to 
over the range of the figure. For these problems Eq. ([10]) gives somewhat lower 



n 



i.i 



search costs than Eq. (J^j. The slow growth in search cost is particularly impressive 
since an unstructured quantum search requires of order 2 n//2 steps, which for n = 100 
is about 10 15 . While this scaling is a polynomial growth in search cost, its cost grows 
a bit faster than that of classical heuristics for these problems. 

For Eq. (|9|), the number of steps giving the smallest cost grows very slowly, ranging 
from 2 to 4 over the range of variables shown in the figure. For Eq. ([II]), the largest 
probability for a solution is after n/2 + 1 steps. However, for the problems with 
n > 200, the small probability of a solution after just the first step is enough to give 
a slightly smaller search cost. 

Further insight into the behavior of this algorithm is given by Fig. ||| which shows 
how the probability to have assignments with different numbers of conflicts varies 
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Figure 2: Expected search cost for maximum constrained 1-SAT vs. n with the two phase 
choices: black and gray are Eq. @ and Eq. (|l0|), respectively. The lines are least squares 
fits to the points. 




conflicts 



Figure 3: Evolution of the probability in assignments with different numbers of conflicts 
for the maximum constrained problem with n = 100 using Eq. (fLC|). The probability for a 
solution reaches about 0.3 after 51 steps, giving a search cost of about 170. 




with each step of Eq. ((^). Specifically, for each step j and each number of conflicts c, 
the figure shows the value of J2 S IV^I 2 where the sum is over all assignments s with 
c conflicts. As described in Appendix p.2j , the amplitudes depend only on the 
number of conflicts in the assignment s. Thus each term in this sum is the same, giving 
(c) l^c 7 ^ 2 where ip^ denotes the amplitude of any of the assignments with c conflicts. 
The initial condition (not shown) has equal probability, 2~ n , in all assignments so the 
corresponding values in this plot would be 0£j2~ n . We see a large probability in states 
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with n/2 — (j — 1) conflicts after step j, illustrating the effectiveness of Eq. (|H]) in 
using neighbor relations to move amplitude. Moreover, the large spike at 50 conflicts 
for step 1 shows the effectiveness of Eq. ( |10a|) in concentrating amplitude in states 
with Abetter = A^start for this problem. In this case, 

(so) l^so I 2 = °- 39 compared to 

its initial value (^) 2 ~ m = °- 08 - 

The relatively low search costs and polynomial scaling are also observed for prob- 
lems with somewhat fewer constraints than the maximum, e.g., m = m max /2, using 
this choice of phases. While this extends the range of good performance for the 
quantum algorithm, such a scaling also gives relatively easy problems for classical 
heuristics as well. 



5.3 Problems with Fewer Constraints 



The ability of this algorithm to concentrate amplitude into solutions for highly con- 
strained problems is a significant improvement over unstructured search methods. 
However, such problems are inherently easy because they can be readily solved by 
classical heuristic methods, for both incremental and repair style searches. By con- 
trast, the difficult search problems, on average, have an intermediate number of con- 
straints [p3[] : not so few that most assignments are solutions, nor so many that any 
incorrect search choices can be quickly pruned. Specifically, for /c-SAT the difficult 
problems have scaling m = 0(n), with the proportionality constant giving the largest 
concentration of hard cases depending on the choice of k. For instance, with k = 3, 
the hard cases are concentrated near [HI m — 4.2n. Thus it is important to examine 
the behavior of the algorithm for problems with fewer constraints. In this section we 
do so using the phase choice of Eq. (0) . 



. 006 
. 005 
. 004 



^0.003 
£ . 002 



. 001 





1 



2 ' 3 ' 4 ' 5 ' 6 ' 1 ' 
conflicts 



Figure 4: Probability in each assignment for a 3-SAT instance with n = 12 after J = 5 
steps. The problem has m = 48, with 4 solutions and -P so in = 0.013. The large gray points 
are the solutions. 
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A larger example of how the algorithm concentrates amplitude into solutions is 
shown in Fig. [|. It shows the values of l^" 7 -*! 2 for each of the 2 n assignments. The 
assignments are grouped by their number of conflicts and, within each group, by the 
integer corresponding to the binary representation of the assignment. For clarity, 
assignments with each number of conflicts are given the same amount of horizontal 
space in the plot, even though there are, e.g., many more assignments with 6 conflicts 
than with or 12. There are 2 12 assignments, so random selection would give a 
probability of about 0.0002 to each assignment. We see that the algorithm results in 
many states with relatively few conflicts, including the solutions, having considerably 
larger probabilities than due to random selection. 

The figure also illustrates the variation in |V4 J ^| 2 among the assignments showing 
that, unlike the extreme problem of the previous section, the amplitudes do not 
depend only on the number of conflicts. Rather the details of which constraints 
apply to each assignment give rise to the variation in values seen here. This variation 
precludes a simple theoretical analysis of the algorithm. 



c 




n 



5 10 15 20 



Figure 5: Expected search cost for random 3-SAT with at least one solution vs. n with 
m/n = 4 on a log plot. For each value of n, 1000 problem instances were used. Error bars 
showing the expected error in the estimate are included but are smaller than the size of the 
plotted points. 

For problems near the transition from over- to underconstrained cases, the scaling 
behavior of this algorithm is shown in Fig. [5]. Specifically, we generated random 
soluble instances of 3-SAT problems with m = 4n, as described in With m = An, 
a large fraction of the randomly generated instances are soluble, so random soluble 
instances are readily generated. The nearly linear behavior on this log-plot indicates 
the seach cost grows exponentially for these problems. Thus this algorithm is not 
particularly effective for the hard problem instances, which are concentrated near the 
transition region of m = 4.2n. 
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Figure 6: Expected search cost for random 3-SAT for n = 10 vs. m/n. Each point 
represents 1000 problem instances, except only 100 random soluble instances are shown at 
m/n equal to 12 and 13. The black line is for random soluble problems while the gray line 
is for problems with prespecified solution. Error bars showing the expected error in the 
estimate are included but, in most cases, are smaller than the size of the plotted points. 

To show how this algorithm is capable of using structure to improve performance, 
we examined the behavior as a function of the number of constraints. To allow sam- 
pling the behavior of problems with many constraints we used random problems with 
prespecified solution to cover the full range of m/n, generated as described in §[| 
For comparison, we also examined the behavior of random soluble problems up to 
m/n = 15. As shown in Fig. [], the search cost eventually decreases as problems 
become highly constrained, though at a larger value of m/n than for both the incre- 
mental quantum search method [^2j and classical methods, whose cost is largest near 
the transition at m = 4.2n. Thus the local quantum search method introduced here 
has only limited effectiveness at using structure to reduce search cost. Nevertheless, 
it is better than unstructured search for highly constrained problems. 

6 Discussion 

The algorithm presented here shows how the structure of search problems can be used 
as the basis of a quantum search algorithm. The algorithm is particularly effective 
for relatively highly constrained problems. It is less effective for problems with an 
intermediate number of constraints. 

The algorithm might be improved in a number of ways. First, the initial motiva- 
tion for the matrix U was to maximally connect assignments to their neighbors. In 
fact, we found that the mapping allowed the algorithm to work best with fewer steps 
than might be expected from moving among neighbors one step at a time. It may be 
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possible to design other mappings that do this even more effectively. 

Another issue is the structure of the types of mappings possible with different 
choices of the phases for the diagonal matrix D. As we saw, the matrix elements 
depend only on the Hamming distance between the assignments when the diagonal 
elements of D depend only on the number of 1-bits. It may be helpful to apply the 
testing operation, and consequent phase adjustments, more frequently as amplitude is 
moved to neighbors than is possible when U has the largest possible mixing between 
neighbors. One such method is to decompose U into a product of matrices, each of 
which introduces a smaller amount of mixing, e.g., using U = WD 1 ' S W, where Z) 1//s is 
a diagonal matrix with Th = 1 if h < n/2 and r = e m ^ s otherwise. Phase adjustments 
can then be introduced between the components of this product. Another method 
uses somewhat smaller values of U\ while retaining real values for U by picking a 
value a < n/2 and using = 1 if h < a and r = — 1 otherwise. The resulting 
mixing matrix has a more diffusive behavior, giving smaller changes with each step 
and hence more opportunity to apply the phase adjustment. Taken to an extreme, 
with a = 0, this recovers the diffusion matrix of the unstructured search |I8|, which 
moves too little amplitude among states at each step to finish rapidly. 

There are also a variety of phase adjustment policies. Those studied here are ef- 
fective for highly constrained problems, but we found that other choices can enhance 
performance in some cases. Furthermore, since we focus on typical or average be- 
havior, other choices that do not improve the average but result in smaller variance 
would also be useful in improving the predictablility of the algorithm's performance. 

As a possible extension to this algorithm, it would be interesting to see whether the 
nonsolution sets with relatively high probability could be useful also. If so, the benefits 
of this algorithm would be greater than indicated by its direct ability to produce 
solution sets. This may also suggest similar algorithms for the related optimization 
problems where the task is to find the best solution according to some metric, not 
just one consistent with the problem constraints. 

There remain a number of important questions. First, how are the results degraded 
by errors and decoherence, the major difficulties for the construction of quantum 
computers [f27| , [40| , |19| , |34|j ? While there has been recent progress in implementation |IJ 



T7\, |39|, quantum approaches to error control 0, |37[ and studies of decoherence 



in the context of factoring |7J it remains to be seen how these problems affect the 
algorithm presented here. However, two aspects of our algorithm may reduce the 
effect of errors. First, the algorithm is based on attempts to randomize phases of 
contributions to nonsolutions while maintaining similar phases for contributions to 
solutions. Since the precise values of the randomized phases are not critical, some 
additional variations due to noise should be tolerable. Second, for problems with n 
variables the algorithm requires only 0(n) steps during which quantum coherence 
must be maintained. That is, coherence need not be maintained between successive 
trials of the algorithm. This contrasts with the 0(2 n//2 ) coherent steps required by 



the unstructured algorithm [18]. Thus even though our algorithm may need to be 
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repeated many times to give a good chance of finding a solution, it has less stringent 
coherence requirements which can simplify its hardware implementation. 

The second remaining question is due to the algorithm's restriction to CSPs with 
two values per variable. This contrasts with the method that constructs solutions 
incrementally |^TJ by operating in a greatly expanded search space. Other CSPs can 
be converted to equivalent problems with two values, but it is often better to operate 
directly with the problem as specified. Thus an important open question is whether 
the neighborhood-based mixing matrix can be effectively generalized to apply directly 
to CSPs with larger domain sizes. 

Third, it would be useful to have a theory for asymptotic behavior of this algorithm 
for large n, even if only approximately in the spirit of mean-field theories of physics. 
This would give a better indication of the scaling behavior than the classical simula- 
tions, necessarily limited to small cases, and may also suggest better phase choices. 
Considering these questions may suggest simple modifications to the quantum map 
to improve its robustness and scaling. 

Finally, there is the general issue of optimally using the information that can 
be readily extracted from search states. Local search methods rely on the number 
of conflicts and the neighborhood relations among states. In the method presented 
here, only a small amount of this information was actually used to determine the 
phase choices. Additional information is available on partial assignments, as used 
with incremental searches, but at the cost of involving a greatly expanded search 
space. Making fuller use of this available information could improve the performance, 
especially in conjunction with a theoretical understanding of the typical structure of 



classes of search problems [23 



A Appendix: Neighborhood Mixing Matrix 

Here we show that the mixing matrix U used in our algorithm is the unitary ma- 
trix depending only on the Hamming distance between states that gives the largest 
possible contribution to neighbors. We do this in two steps. We first show that matri- 
ces whose elements depend only on Hamming distance can be expressed in the form 
WDW where W is given by Eq. (|7|) and D is diagonal with elements depending only 
on the number of 1-bits in the assignments. In the second step, we show that our 
choice of values for D gives the largest possible mapping to neighbors among unitary 
matrices of this form. 

A.l Matrix Decomposition 

For n variables, the matrices have dimensions 2™ x 2 n and elements determined by 
a set of values uq, . . . , u n , i.e., U rs = Ud( r , s )- For example, with n = 2 and the states 
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binary order, i.e., 00, 01, 10, and 11, the matrix has the general form 



U 



( U Ui Ui u 2 \ 

Ui U U 2 Ui 

Ui U 2 U Ui 

\u 2 Ui Ui U J 



(11) 



These matrices have a simple recursive decomposition when the states are in 
binary order, i.e., ordered by the value of the integer with corresponding binary 
representation, namely 

where the 2 n_1 x 2 n_1 matrices Uq and U\ also have elements depending only on the 
Hamming distance between assignments, but considering only the first n— 1 variables. 



U 



(12) 



An example of this decomposition can be seen in the structure of Eq. (11). 

This decomposition gives a particularly simple expression for these matrices. 
Specifically, define the matrix W as 



ArAs\ 



(13) 



For example, when n = 2, 



(I 1 1 


1 \ 


1 -1 1 


-1 


1 1 -1 


-1 


[l -1 -1 


1 / 



w 

This matrix also has a recursive decomposition 

W=[ W> -W> 



where W ' is the same matrix but defined on assignments to n 
Now consider the product WUW. For n = 1, this gives 



(14) 



(15) 



1 variables. 



2(w + u x ) 






2(u - ux) 



(16) 



which is a diagonal matrix. Now suppose that WUW is always diagonal for any 
matrix with U rs depending only on d(r, s) up to n — 1 variables. Using the recursive 
decompositions of Eq. (|I~2"D and (|T5|) then gives, for a matrix with n variables 



WUW 



2W'{U + Ui)W 






2W'{U - UJW' 



(17) 



Each of the submatrices are diagonal because they involve products of the form WUW 
for n — 1 variables. Thus, by induction, WUW is a diagonal matrix for all n, which 
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we denote by ND, where N = 2". Note that Eq. (0) is W = W/y/N. The matrix W 



is symmetric and unitary ||18|| . In particular, W is its own inverse: WW = I. Thus 
from WUW = ND we obtain U = WDW. 

Furthermore, the diagonal elements of D depend only on the number of 1-bits in 
the assignments, i.e., D rr = r\ r \ for some set of values To, ... , r n . To see this, we have 

(WUW) rr = £^(^)(- 1 ) |rAX|+bAr| (18) 
xy 

\r\ 



E^Et-^E' 1 (is) 



d 2=0 xy 

where the innermost sum of the second line is over all assignments x and y for which 
d(x, y) = d and |r A x\ + \y A r| = z. For this inner sum, consider those assignments 
x with exactly 77 1-bits in common with r, i.e., |r A x\ = 77. There are ^^2 ri ~' r ' 
such choices for x. Then the y sum counts the number of assignments y such that 
d(x, y) = d and \y A r\ = z — rj. These assignments are characterized by a, the number 
of 1-bits they share with both r and x. The number of such y is readily seen to be 




Summing over a then gives the total number of assignments y, which we denote as 
v(d,z,r], \r\). Combining these observations gives 



(WUW) rr = E^E(-!) Z E ( )2 n - l Md,z }V} |r|) (20) 
d 2 n V V 



which depends only on the number of 1-bits in the assignment r. 

In summary, any matrix whose values depend only on the Hamming distance 
between assignments can be written in the form U = WDW where D is a diagonal 
matrix with elements depending only on the number of 1-bits in the assignments. 
That is, the values for u , . . . ,u n are uniquely determined by the diagonal values of 
D, which in turn are fully specified by the values r , . . . , r n for assignments with each 
possible number of 1-bits. 

Since W is unitary, U will be unitary if and only if D is. For diagonal matrices, 
the unitary condition is equivalent to the requirement that each diagonal element is 
a phase, i.e., a complex number whose magnitude equals one. 

A. 2 Maximum Mapping to Neighbors 

To determine the choices for r h that give the largest possible value for ui, we have 

1 n 

U rs = ^T, T hS h (r,s) (21) 
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where 

S h (r,s)= £ (-l) lrMl+lsMl (22) 

t,\t\=h 

with the sum over all assignments t with h 1-bits. 

A given 1-bit at position e of t contributes 0, 1 or 2 to |r A t\ + \s A t\ when positions 
e of r and s are both 0, have exactly a single 1-bit, or are both 1, respectively. Thus 
(_l)|rAt[+|*At| e q Ua i s (— l) 2 where z is the number of 1-bits in t that are in exactly one 
of r and s. There are (|r| — |r A s|) + (|s| — |r A s\) positions from which such bits of t 
can be selected, and by Eq. (D) this is just d(r, s). Thus the number of assignments t 
with h 1-bits and z of these bits in exactly one of r and s is given by m (jjlf) where 
d = c?(r, s). Thus S^r, s) = where 

so that [/ rs = w rf ( rjS ) with u d = E^S^ /N. 

To select the values of that maximize u±, note that Sqi = 1 and 5*^™^ = 
l^™^ 1 ^ — Uli)- Thus is positive for n > 2h and negative for n < 2h, and iti is 

maximized by selecting to be 1 for h < n/2 and —1 for h > n/2. If n is even, 

is zero for h = n/2 so the choice of r„/ 2 does not affect the value of u\. In this case, we 

take T n / 2 = 1. These choices give U\ = ^(i^T^j) which scales as ^2/ (nn) as n — > oo. 
Note this is much larger than the off-diagonal matrix elements in the diffusion matrix 



used in the unstructured search algorithm IIS], which are 0(1/N) =0(2 



This maximum neighbor mixing matrix is identical to that used with an incre- 



mental search method [25| but with a different interpretation of the search states. 
If instead we take tq = — 1 and the remaining values equal to one U becomes the 
diffusion matrix of the unstructured search method [|18j. 



B Appendix: Classical Simulation 

A direct classical simulation of the quantum algorithm results in an exponential in- 
crease in the cost, limiting the empirical evaluation of the algorithm to small cases. 
Nevertheless, the simple structure of the mixing matrix U can be used to reduce this 
cost penalty. It can be further reduced for the special case of soluble problems with 
the maximum possible number of constraints. 

B.l General Problems 

As a practical matter, it is helpful if a quantum search method can be evaluated 
effectively on existing classical computers. Unfortunately, the exponential slowdown 
and growth in memory required for such a simulation severely limits the size of feasible 
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problems. For example, Eq. ([5]) is a matrix multiplication of a vector of size 2™ so a 
direct evaluation requires 0(2 2n ) multiplications. 

The cost of the classical simulation can be reduced substantially by exploiting the 
map's simple structure. Specifically, the product F^x can be computed recursively 
using Eq. (|T^) giving 

W * = [wW-W'x?)) (24) 

where x^ 1 ) and x'- 2 -' denote, respectively, the first and second halves of the vector x. 
Thus the cost to compute Wit is C(n) = 2C{n — 1) +0(2 n ) resulting in an overall cost 
of order n2 n for this product as well as the full mapping step ?7x = (l/N)WDW~x.. 
While still exponential, this improves substantially on the cost for the direct eval- 
uation on classical machines. An open question is whether other techniques could 
give approximate values for the behavior without an exponential cost on classical 
machines ||. 



B.2 Maximally Constrained Problems 

In general, the cost of a classical simulation of the quantum algorithm grows exponen- 
tially with the number of variables. However, the simulation cost is greatly reduced 
in some special cases where the amplitudes have a regular structure. One such case 
is provided by problems with the maximum possible number of constraints that still 
allows for a solution. 

For such problems, the amplitudes depend only on the number of conflicts in 
each assignment, not their particular location. This observation allows for a smaller 
representation of the search state and hence an empirical study of larger cases. To 
see this, consider a 1-SAT problem with n variables. The most constrained, but 
still soluble, case will have n conflicts, each involving a distinct variable. Thus each 
variable has an associated "bad bit" value that causes a conflict. The number of 
conflicts in a given assignment s is then just equal to the number of bad bits it 
contains. An assignment with c conflicts will have c neighbors with c — 1 conflicts and 
n — c neighbors with c + 1 conflicts. Thus the neighborhood structure of the problem 
is uniquely determined by the number of conflicts in the assignments. 

Suppose ij}^ depends only on the number of conflicts in the assignment s. We 
need to show that after a single step, Eq. (H) gives values for that only depend 

on the number of conflicts in r. First note that the phase choice p s uses only the 
neighborhood structure of s, which for these problems depends only on the number 
of conflicts in s. Thus p s ips depends only on c, the number of conflicts in s, and 
can be denoted by <p c . Since the mixing matrix U rs depends only on the Hamming 
distance between the assignments r and s, Eq. (|) becomes 

^ +1) = E^E^E'i (25) 

d c s 
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where the inner sum is over assignments s with c conflicts and Hamming distance d 
from the assignment r. Suppose r has b conflicts. Of the c "bad bits" in s, (c+b—d)/2 
are in common with those of r, and the remaining (c — b + d)/2 do not appear in r. 
The number of ways to construct such assignments s is then 




Thus the result of this mapping step depends only on the number of conflicts in 
the assignment r, preserving the simple structure of the search state. Thus we can 
represent the entire search state by simply keeping track of the amplitudes associated 
with each possible number of conflicts, i.e., 

= E Vg~Pc4><P (26) 

c 

where b and c range from to n and 



d 



Xc + b-d)/2) \(c-b + d)/2 



Because the search state depends only on the number of conflicts and not their 
specific associated variables or values, we can, without loss of generality, choose the 
conflicts so that all the "bad bits" have the value 1, i.e., the unique solution of the 
problem has all variables equal to 0. With this choice, we can then make use of the 
decomposition of the mixing matrix. That is, V max = jy max £) max yj/ max where _D max 
is a diagonal matrix with -D^ ax equal to 1 for b < n/2 and —1 otherwise, and 

«=-=^?(-^QK)=S (28) 



from Eq. (f23|). Here the binomials in the sum count, for an assignment r with b 1-bits, 
the number of assignments s with c 1-bits that have |r A s\ — z. This decomposition 
is possible for this particular choice of the unique solution because the number of 
1-bits correspond to the number of conflicts. 

In this way, classical evaluation of Eq. ([5]) reduces to multiplication by matrices 
of size (n + 1) x (n + 1), with cost of order n 2 . 

Finally, 1-SAT problems with fewer than the maximum number of constraints also 
have the property that the amplitudes depend only on the number of conflicts in each 
assignment. So this compact representation could be used to study other cases of 
1-SAT with many more variables than is feasible for a direct classical simulation. 
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